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After a brief overview of the realization of chiral symmetry in hot and/or 
dense medium, the implications of the existence of a light scalar-isoscalar 
meson (which we call a throughout this report) are discussed in the vacuum 
and in the medium. Special emphasis is put on its relation to the fluctuation 
of the chiral order parameter (qq). In the vacuum, the c-meson is an elusive 
resonance corresponding to a pole deep in the second Riemann sheet of the ir- 
7r scattering matrix in the I=J= 0 channel. This is because the the amplitude 
fluctuation of (qq) (corresponding to ur-c) and the phase fluctuation of (qq) (2 
pion states) mix strongly in the vacuum. As the temperature and/or density of 
the system increase, however, there arises a softening of this complex pole due 
to the partial restoration of chiral symmetry. Such soft modes play a key role 
to understand the QCD phase structure. We demonstrate that even a slight 
softening of the a mode could induce strong spectral enhancement and strong 
7T-7T correlation near the 2m n threshold in the I=J= 0 channel, which thereby 
can be a signature of the partial restoration of chiral symmetry. Such spectral 
enhancement may be seen not only in hot matter created by the relativistic 
heavy ion collisions but also in cold matter (heavy nuclei) probed by photons 
and hadrons. Relevance of the partial restoration of chiral symmetry in heavy 
nuclei with the recent data by the CHAOS and CB collaborations as well as 
with on-going and future experiments are also discussed. 
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I. INTRODUCTION 


One of the most intriguing phenomena in quantum chromodynamics (QCD) is the dy¬ 
namical breaking of chiral symmetry (DBCS). This explains the existence of the pion and 
governs most of the low energy phenomena in hadron physics. DBCS is associated with the 
condensation of quark - anti-quark pairs in the QCD vacuum, (qq), which is analogous to 
the condensation of Cooper pairs in the theory of superconductivity |T|. As the temperature 
(T) and/or the baryon density (%) increase, the QCD vacuum undergoes a phase transition 
to the chirally symmetric phase where (qq) vanishes. Studying the mechanism of DBCS, 
exploring the phase structure of the QCD ground state, and finding possible signatures of 
the QCD phase transition are the central issues in modern hadron physics 0. 

Various phases of QCD are characterized not only by T and ub, but also by the current 
masses of light quarks ( m U} d , s )• In the real word, the masses of u and d quarks are much 
lighter than the s quark 

m m 

— = 18.9 ±0.8, — = 0.553 ±0.043, (1) 

m d rn d 

where m s (2GeV)(~ 100 MeV) is comparable to A q C d (the QCD scale parameter) and also 
to the critical temperature (density 1 / 3 ) of the quark-hadron transition. Therefore, studying 
the phase structure in the m s -T-riB space by neglecting m u> d as a first approximation would 
give us a good insight into the real world. 

Shown in Fig.l is a possible phase structure projected onto the m s -T plane with ns— 0 
(the left panel) and onto m s -riB plane with T=0 (the right panel). m s =0 (m s = oo) corre¬ 
sponds to the limit of SUl(3 ) <E> SUr( 3) (SUl( 2) <8> SUr(2 )) chiral symmetry. 



FIG. 1. Possible phases in the m s -T plane at ur =0 (left panel) and the m s -nB plane at T=0. 
m U)d = 0 is assumed. Relevant degrees of freedom in each phase are also shown, no = 0.17fnU 3 is 
the normal nuclear matter density. 
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The chiral phase transition for T ^ 0 and Ub = 0 has been extensively studied on the 
lattice for recent years with (uu)(= (dd)) as the chiral order parameter Q. The analysis 
based on the renormalization group and the universality hypothesis indicates that the chiral 
transition at finite T is of second (first) order for large (small) m s This has been con¬ 
firmed by the lattice QCD simulations with dynamical quarks. The critical temperature for 
massless 2 and 3 flavors has been also studied on the lattice and the results are summarized 
as 


T pc (N f = 3) = (154 ± 8)MeV 
T pc (N f = 2) = (173 ± 8)MeV 
T pc (N f = 2) = (171 ± 4)MeV 


(staggered fermion), 
(staggered fermion), 
(Wilson fermion), 


( 2 ) 


where T pc is a pseudo-critical temperature extracted from the chiral susceptibility. 

The second order line and the first order line in Fig .[I] (left panel) is expected to meet 
at some intermediate value of m s at the so-called tricritical point || . However, the precise 
location of this tricritical point in m s -T plane is not identified yet. If we have finite but small 
u, d quark masses, the second order line turns into a smooth crossover, and the tricritical 
point becomes an end point of the first order line. Therefore, in the real world where two 
light quarks + one medium-heavy quark exist, the chiral transition at T ^ 0 is either the 
first order or the crossover. 

The chiral phase transition at finite baryon density is not well understood. One of the 
reasons is the first principle lattice QCD simulation for finite chemical potential (//) has 
a severe sign problem originating from the complex fermionic determinant [|4|. Analyses 
based on simple models such as the Nambu-Jona-Lasinio (NJL) model show that the chiral 
transition may be first order for ^ 0 and T -C T c 0. [] If this is the case, the quark- 
hadron mixed phase is expected in the m s -nB plane as schematically shown in the right panel 
of Fig.l. The system at finite baryon density, which is intrinsically quantum, is obviously 
richer in physics than that at finite T and has various phases such as the baryon superfluidity, 
meson condensations, and the color superconductivity (see the reviews, 0). 

For sufficiently large m s with m u ^ = 0, there will be also a tricritical point in the T-/i 
plane. For small but finite m Ut d, the tricritical point again becomes an end point of the 
first order line. Very recently, a new attempt on the lattice has been made to locate this 
end point by looking for the Lee-Yang zero of the partition function numerically, where the 
2-dimensional reweighting method is employed to obtain the partition function at finite T 
and // [|I0| . It is found that T en d=(160±3.5) MeV and /i en d=(725±35) MeV for (2+l)-flavors 


using the data on 4 3 x 4, 6 3 x 4 and 8 3 x 4 lattices. Although obtaining a clear signal of 
Lee-Yang zeros may require much larger lattice size, this approach provides a new way of 
analyzing the physics near the critical line in the T-/i plane. 


3 The possible vector coupling between quarks in the NJL model may alter the order of the chiral 
transition at finite density with low temperatures. [|J 
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II. THE LIGHT SCALAR-ISOSCALAR MESON IN QCD 


The order parameter (qq) of the chiral transition is determined as the value at which the 
effective potential (free energy) V(cr) takes the minimum. The light scalar-isoscalar meson, 
which we call o in the following, is a particle representing the amplitude fluctuation of the 
order parameter around the minimum of V(a ). In this sense, the a meson is analogous to 
the Higgs particle H in the Glashow-Salam-Weinberg model. However, the properties of 
cr and H or, more precisely, the physical modes in the respective channels, could be quite 
different because of the following reason: The Nambu-Goldstone (NG) boson, which is a 
phase fluctuation of the order parameter, appears as the pion in QCD, while the NG boson 
in the electro-weak theory is absorbed into the gauge bosons. Therefore, a can have a strong 
width decaying into two pions, while such process for H does not exist. In other words, the 
physical a state in QCD can be at most a broad resonance and is represented by a linear 
superposition of the ur-cr (genuine amplitude-fluctuation of the order parameter) and the 
two-pion state. 


A. Chiral symmetry, analyticity, crossing symmetry 
and complex <r-pole in n-n scattering 


As mentioned above, a tricky point with the a meson is that it couples to two pions to 

firm 


acquire a large width T^ ~ Re m a and make itself elusive 


Nevertheless, recent careful 


phase shift analyses of the tt-tt scattering in the I—J —0 channel have come to claim a cr pole 
in the complex energy plane with the real part Re m a = 500-800 MeV and the imaginary 
part Im m a ~ 500MeV 


12.11.13 


One of the crucial points to deduce the poles in the -k-tt scattering matrix is to construct 
the invariant amplitude so that it satisfies the chiral symmetry constraints, the analyticity, 
the unitarity and the (approximate) crossing symmetry. For this purpose, some reliable 
resummation method such as the N/D method should be employed. This point has been 
emphasized e.g. in [113] and in |T5|,[T6[|. In the construction of 


the ur-a strongly coupled 
to 27 t is favored to have a good overall agreement with the phase shifts in I=J= 0, I—J =1 
and (1=2, J— 0) phase shifts. In |T5yi6|] , the unitalization of the strong 2n correlation in 
the I=J =0 channel is enough to reproduce the phase shifts. Although it is not clear at the 
moment whether one really needs to introduce the ur-a or not to reproduce the phase shifts, 
the existence of the light cr pole with a large imaginary part in the 2nd Riemann sheet is 
the common feature in all modern analyses. This is seen in Fig.|2] where a recent summary 
(taken from |17]]) of the complex a pole is shown. 
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FIG. 2. The cr-pole of the 7r-vr scattering matrix in the complex energy plane in various analyses. 
This figure is taken from [|T7|]. 


B. Hadron phenomenology and the a meson 


If the location of the a pole in the 2nd Riemann sheet has Re m a ~ 500 MeV, many 
experimental facts which otherwise are mysterious can be nicely accounted for in a simple 
way. Some examples in hadron phenomenology where the scalar-isoscalar fluctuations are 


relevant include the following. More extensive discussions can be see in |]TT 


The scalar-isoscalar resonance with Re m a ~ uik may have a sizable contribution to 
the enhancement of the A/ = 1/2 process in K— > mr decay ||TB . 


The state-independent attraction, which is responsible for the nuclear binding, orig¬ 
inates from the two-pion exchange potential. The two-pion exchange includes the 
ladder , the cross and the rescattering diagrams with the A(1232) being incorporated 
in the intermediate states. The a plays a role through the rescattering contribution as 
well as through the possible direct a-N coupling. The former should be constructed 
consistently with the tt-tt phase shift in the I=J =0 channel. For the latter, systematic 
analysis of N-N, hyperon-jV and hyperon-hyperon interactions may be useful Jl9j . 


The tt-N sigma term and the generalized scalar-charge of baryons can be defined 
as Sjrjv = m(N\uu + dd\N) and (B\qq\B) (q = u,d,s,...), respectively. Effective 
charges are usually enhanced (suppressed) due to collective excitations generated by 
the attractive (repulsive) forces. The enhancement is caused by the polarization of 
the vacuum in the I=J =0 channel and can indeed account for the empirical value of 
— 45 MeV without introducing a large OZI violation in the nucleon [|20| . 


There are some effective theories of QCD which predict the light ur-c such as the 
Schwinger-Dyson approach with the rainbow-ladder approximation pi], the Nambu- 


Jona-Lasinio model p2], and an approach based on the mended symmetry p3 . 
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The above list implies that the correlation in the scalar-isoscalar channel is indeed im¬ 
portant in the hadron dynamics. This is in a sense natural because the dynamics which is 
responsible for the correlations in this channel is nothing but the one which drives the chiral 
symmetry breaking. 

A fundamental question is, then, whether the complex cr-pole observed in the 7r-7r phase 
shift can be interpreted as a quantum fluctuation of the order parameter associated with 
the dynamical breaking of chiral symmetry or its remnant. In condensed matter physics, 
such a question can be answered by changing the ground state properties of the system with 
temperature, pressure, doping etc as external controlling parameters J23J. If the frequency 
of some collective mode has characteristic change near the critical point, it provides us with 
a clear signature of a direct link between the mode and the phase transition. The life is, 
unfortunately, not that easy in QCD from the experimental point of view. Nevertheless the 
heavy nuclei provide us with finite baryon-density environment which allows us to study the 
partial restoration of chiral symmetry and associated change of the meson properties. Heavy- 
ion collisions also provide us with not only high baryon density but also high temperature 
although the system stays in such environment only in a transient time. 

Some years ago, one of the present authors p5[ proposed several nuclear experiments 


including one using electro-magnetic probes for producing the a meson in nuclei to see 
a clearer evidence of the a pole and to explore possible restoration of chiral symmetry 
in nuclear medium. When a particular meson is put in a nucleus, it may dissociate into 
complicated excitations to loose its identity in the medium; for example, a 2i r, p-h, n+p- 
h, A -h, n+A-h, ■ ■ ■. Therefor it is not appropriate to talk about the “mass” and “width” of 
the particle, but one should study the spectral function which contains all those information. 


III. CHIRAL RESTORATION AND SOFT MODE AT FINITE T AND DENSITY 

The basic observation on which the whole discussions in this report are based is that 
the dynamical breaking of chiral symmetry is a phase transition of the QCD vacuum with 
an order parameter ( qq ): Then there may exist collective excitations corresponding to the 
quantum fluctuations of the order parameter. The phase (amplitude) fluctuation of the 
order parameter corresponds to the pion (ur-cr). As we have mentioned before, they couple 
strongly and mix through the process a 2 tt. 

Now let us first show how the magnitude of the condensate (qq) decreases at finite T 
and/or density Ub in an analytic way. 


A. Chiral condensate at low T 


The condensate at T ^ 0 (n B = 0) can be written as 


(qq) 


-Tr 


qq e 


Hqcd/T 


df(T) 

dm q 


dP(T) 

dm q 


( 3 ) 


where Z is the QCD partition function, / is the free energy density and P is the the total 
pressure of the system including vacuum pressure and the thermal pressure. This is because 
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the current quark masses enters Hqcd only in the form J2i m qQQ] hence this formula is valid 
for each flavor (q = u, d, s, c, b , t). 

At high T A qcd, one can calculate P(T ) using the high T QCD perturbation theory 
and (qq) is shown to be exactly zero for m q = 0. This is because there is no vertex which 
flips chirality in the QCD Lagrangian except for the mass term. On the other hand, there 
is a dynamical breaking of chiral symmetry at low T and (qq) does not necessarily vanish. 
For T low enough, the system may be approximated by a dilute and weakly interacting gas 
of pions. Its pressure can be calculated using the chiral perturbation theory |[26|| : 


Z = Tr 


,-Hqco/T 


exp 


r-l/T 


dr / d?x C e ff(U) 


( 4 ) 


where C e ff(U) is expanded as C^ + C (6 ) + • • ■ in powers of <9/4 irf n and rrin/^fn- 

Therefore, the quark mass dependence of P(T) for Nf = 2 at low T in (|3j) can be extracted 
from its m^ dependence through the Gell-Mann-Oakes-Renner relation 




= -2 + 
J 7r 


( 5 ) 


where m = (m u + rrid)/2 , and (qq )o and f n are the vacuum condensate and the pion decay 
constant in the chiral limit, respectively. 

Writing P(T ) = P p i 0 n(D) + P vac , and using the above relation, one thus finds the quark- 


condensate in the chiral limit for Nf = 2 as 

T 2 1 (T 2 


(QQ) = [+ 1 dPpion (T) 


(qq) 


fn dml 


= 1 - 


m-n^O 


m e v s fi 


f (5) h (T) +o(r8 >' <6) 


where A q (= 470 ±110 MeV) is a parameter extracted from the experimental tt-ti scattering 
length in the isoscalar D -wave channel. The T 2 corrections in the above expansion originate 
solely from the free pion-gas contribution. Therefore, it can be alternatively obtained by 
applying the soft pion theorem to the free pion gas. The negative contribution of the 0(T 2 ) 
term is due to the fact that the pion has positive scalar-charge (n(p)\qq\'K(p)), which can 
be seen from either by the soft-pion theorem or by the Feynman-Hellman theorem. On the 
other hand, the interaction among pions are relevant for the terms 0(T n - v ). 

The above result of the low T expansion suggests that the chiral condensate decreases 
uniformly as T increases and the symmetry restoration takes place somewhere around 100- 
200 MeV. However, the behavior of qq near the critical point is out of reach of the low T 
expansion, because the system is not any more the dilute gas of pions near the critical point. 
In such a region, non-perturbative methods such as the the lattice QCD simulations are 
necessary to make a reliable estimate of the chiral condensate. 


B. Chiral condensate at finite baryon density 

(qq) at finite baryon density obeys an exact formula in QCD directly obtained from the 
Feynman-Hellman theorem 
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( 7 ) 


which is valid for any eigenstates |T) of the QCD Hamiltonian Hqcd(= J cPx'Hqcd ) with 
the normalization (TIT) = 1. For cold nuclear matter T) = |nm), the energy density is 
given as 


(nm | Hqcd | nm) ^vac + n b[M n + E hin ], 


( 8 ) 


where £ vac is the vacuum energy density, Mat is the rest mass of the nucleon and £bin. is the 
binding energy of nuclear matter per particle. Applying this and the similar formula for the 
nucleon to (0), one arrives at [|7|: 


(uu + dd) 
(uu + dd) o 



A di E \» u . (ti ) , 


( 9 ) 


where E^at = 45 ± 10 MeV is the pion-nucleon sigma term. The density-expansion of the 
r.h.s. in the lowest order gives a reduction of almost 35 % of ( qq } already at the nuclear 
matter density no = 0.17fm -3 : The physical origin of this reduction is similar with the pion 
gas case discussed above; the nucleon has the positive scalar charge (N\qq\N) > 0, hence 
the nucleon gas gives the positive contribution to the quark condensate. The above estimate 
suggets that heavy nuclei can be a good laboratory to study the partial restoration of chiral 
symmetry in nuclear medium; 20 to 30% reduction of the chiral condensate might be large 
enough to manifest typical phenomena of partial restoration of chiral symmetry. We shall 
argue that this is actually the case. 


C. Hadronic correlations and symmetry restoration 

The wisdom of many-body theory tells us that the structure of the ground state is closely 
related with and reflected in the excitation spectra of the system. In the case of chiral 
symmetry restoration, a model-independent statement on this connection goes as follows. 
Consider the hadronic correlations in a same chiral multiplet, such as S (scalar)— P a (pseudo¬ 
scalar) and V ^ (vector)— A a (axial-vector). Then they must be degenerate when (qq) —» 0, 
namely, 

(S(x)S(y)) - (P a (x)P a (y)), (A«(x)A b u (y)) - (V^x)V b (y)). (10) 

To see such chiral degeneracy by lattice QCD simulations, it is useful to define the 
thermal susceptibilities for hadronic operators O defined in the Euclidean space, 

X H = / o dr j d3x (0 f (r,x)£>(0,0)). (11) 

Shown in Fig.[3j is the lattice QCD simulation of ^l/x H with Nf = 2 dynamical quarks 
0 . One can see the degeneracy between the cr-channel (7=0, J p =0 + ) and the 7r-channel 
(1=1, J p = 0~) above the critical point. Also, yjl/x H f° r a smaller than that for S at low 
T, which indicates that a could be light. There remains a splitting between the 5-channel 
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(1=1, J p = 0 + ) and the 7r-channel even above the critical point, which reflects the breaking 
of Ua{ 1) symmetry in the 2-point correlation function for Nf = 2; in other words, the 
[/a( 1) anomaly survives the chiral restoration. In general, the effect of the 1 / 4 ( 1 ) symmetry 
breaking appears for N(> Nf)- point function for iVy-quarks as shown in 



FIG. 3. Thermal susceptibilities in three different channels ( 7 r, < 7 , and 5) for two-flavor QCD on 
the 8 3 x 4 lattice with rn u j a = 0.02 ||]. The vertical (horizontal) axis denotes ^ 1 /x H (the lattice 
coupling (5 = 6 /g 2 ). 

To make a direct connection of the chiral restoration with the experimental observables, 
the Euclidean correlations such as ( |TT| ) are not useful enough; one needs to extract the 
information of the spectral distribution as we mentioned before. Recently, a new approach 
has been proposed to extract hadronic spectral functions (SPFs) from lattice QCD data 

0 


by using the maximum entropy method (MEM) 


MEM has been successfully applied 


for similar problems in image reconstruction in crystallography and astrophysics, and in 
quantum Monte Carlo simulations in condensed matter physics p0| . 

The Euclidean correlation function D(t) of an operator 0(r,x) and its spectral decom¬ 
position at zero three-momentum read 


D(r > 0) = 


r, x)O(0, 0))d 3 x = / K(r,u)A(uj)du, 
Jo 


( 12 ) 


where w is a real frequency, and A(u) is the spectral function (or sometimes called the 
image), which is positive semi-definite. K(t,lj) is a known integral kernel (it reduces to 
Laplace kernel e~ TUJ at T=0.) Monte Carlo simulation provides D{ji ) on the discrete set of 
temporal points 0 < r*/a < iV T . From this data with statistical noise, we need to reconstruct 
the spectral function A{u) with a continuous variable aThis is a typical ill-posed problem, 
where the number of data is much smaller than the number of degrees of freedom to be 
reconstructed. MEM is a method to circumvent this difficulty through Bayesian statistical 
inference of the most probable image together with its reliability. There are three important 
aspects of MEM: (i) it does not require a priori assumptions or parameterizations of SPFs, 
(ii) for given data, a unique solution is obtained if it exists, and (iii) the statistical significance 
of the solution can be quantitatively analyzed. 

I 11 Fig.(|. shown is an example of the spectral function in the pion and the rho meson 
channels at T = 0 extracted from the quenched lattice QCD data fl29|. Correct resonance 
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peaks at low energies and the continuum structure at high energies are well reproduced. 
The long-standing problem of in-medium spectral functions of vector mesons (p, tv, (j), <//'*/’, 
T, • • •, etc) and scalar/pseudo-scalar mesons (a, 7r, • • •, etc) can be studied using MEM 
combined with finite T lattice simulations. 


60 

50 

40 
3, 30 

Q? 

20 

10 

0 

0 1 2 3 4 5 6 7 

co (GeV) 










(a) PS Channel 
lattice kernel 
k = 0.1557 






























— 













(b) V Channel 
lattice kernel 
k = 0.1557 

























0 1 2 3 4 5 6 7 

03 (GeV) 


FIG. 4. Spectral functions p ou t( w) = A(lo)/lo 2 obtained by MEM using the quenched lattice 
QCD data. The lattice size is 20 3 x24 with a = 0.0847 fm. 12 data points in the temporal direction 
(1 5; Til a < 12) are used for the MEM analysis, (a) is for the pion channel and (b) is for the 
rho-meson channel. The figures are taken from [p9j]. 


IV. SCALAR CORRELATION IN THE MEDIUM 
A. Spectral enhancement in hot matter 

The fluctuation of the order parameter becomes large as the system approaches the 
critical point of a second order or to the end point of the first order phase transition. This is 
easily seen for the second order case by looking at the Landau free-energy of a double-well 
type written in terms of the order parameter cr rs./ (qq), 

V{°) = _ | (j2 + J (t4 ’ ( 13 ) 

where a(oc T c — T) changes sign at the critical point, while b remains positive. The minimum 
value and the curvature at the minimum of the effective potential <7o read 

= a , 

cr =(JQ 

respectively. As a becomes small, not only the order parameter uq but also the curvature 
decrease. This is nothing but the softening of the oscillational mode of the order parameter 
associated with the second order phase transition. 



fa 1 d 2 V(c r) 

ff0= V6’ 2 da 2 


11 





































1. Soft mode in the Wigner phase 


Historically, the dynamical fluctuation of the order parameter near the chiral critical 
point was first considered in the Wigner phase above T c by using the Nambu-Jona-Lasinio 
model |3TH33ll. It was shown that such collective oscillation becomes soft and narrow as 
the system approaches to the critical point. Since this is the collective mode in the chirally 
symmetric phase, it may be called “para-pion” or “para-sigma” in analogy to the para- 
magnon in the electron gas and liquid 3 He [B4| and para-superconductivity in low dimensional 


superconductors [[3^]. In nuclear physics, the soft mode for the superconductivity is called 


the pairing vibration, and the 2 + phonon is the soft mode for the quadrupole deformation 
of nuclei. |3(J The soft mode for the pion condensation |J should be the longitudinal spin- 
dependent isovector modes with finite momentum |37].|9[]. In Fig.^j. shown is the spectral 
function of the para-pion defined by 


S ab (LU. p; T) oc Im / <fx e* px (R gr a g(x)gT b g(0)), 


(15) 


where T a = (1 ,z7 5 t) and R denotes the retarded product |32 . 



FIG. 5. The spectral function of degenerate “para-pion” and “para-sigma” at finite T in the 
Wigner phase calculated by the 2-flavor Nambu-Jona-Lasinio model in the chiral limit. Criti¬ 
cal temperature of the chiral transition is 164 MeV and the spectral enhancement is seen as T 
approaches to the critical point. Taken from . 


2. Soft mode in the Nambu-Goldstone phase 


In the chirally broken phase just below T C1 there arises two types of soft modes, namely 
the phase and amplitude fluctuations of the chiral order parameter ( qq) 


331. Because of 


the Nambu-Goldstone theorem, the pion stays massless in the chiral limit below T c , while 
one can expect sizable softening (the red-shift) of a. In the real world, the situation is not 
that simple, since a has a large width from the decay a —> 2n. Nevertheless, cr may appear 
as a narrow resonance near T c because of the phase space suppression of the above decay 
channel, although it is at best a broad resonance in the free space. This idea has been put 
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forward in |38|] for the first time and various phenomenological applications in relation to 
the relativistic heavy ion collisions have been discussed later |39| . 

Since the character of a changes from a very broad resonance to a sharp resonance, it is 
not at all enough to just talk about the mass and width but is necessary to study the spectral 
function. By using the linear cr-model, the spectral functions in n and a channels have been 
studied at finite T and it was found that there is a significant spectral enhancement just 
above the 2m n threshold as T increases m- 

To describe the general features of this spectral enhancement, let us consider the prop¬ 
agator of the scalar-isoscalar cr-meson at rest in the medium: 


D„ {ui)=io 


(16) 


where m a is the mass of the a in the tree-level, and T, a (uj]T) is the loop corrections in 
the vacuum as well as in the hot medium. The corresponding spectral function is given by 
p c (uj) = — ^ImZibrju;). Near the 2m n threshold, the imaginary part in the one-loop order 
reads 


ImEo- oc 9(uj 


2 m n ) 



(17) 


When chiral symmetry is being restored, m* (“effective mass” of a defined as a zero of 
the real part of the propagator ReD“ 1 (a; = m*) = 0) approaches to m n . Therefore, there 
exists a temperature at which ReD'^o; = 2 m n ) vanishes even before the complete cr-n 
degeneracy is realized; ReD~ 1 (o; = 2/71^) = [a; 2 — m 2 — ReE cr ] a;= 2 mw = 0. At this point, the 
spectral function can be solely represented by the imaginary part of the self-energy; 


p a (u ~ 2m n ) 


1 

7T IinEo- 


9(u — 2m. K ) 



(18) 


which shows an enhancement of the spectral function at the 2tt threshold. One should note 
that this enhancement is due to the partial restoration of chiral symmetry and hence generic 

4 |. 

To make the argument more quantitative, let us evaluate p a { uj) in the 0(4) linear cr- 
model: 

1 2 \ 

C M = -tr [dMdM^ - p 2 MM ] - — (MM 1 ) 2 - h(M + M% (19) 


where tr is for the flavor index and M — a + ir ■ tt. Although the model has only a limited 
number of parameters and is not a precise low energy representation of QCD at zero T, it 
can describe the pion dynamics qualitatively well up to lGeV HTJ. It also provides a good 


description of the static critical phenomena of QCD owing to the universality argument 
01- The coupling constants /i 2 , A and h have been determined in the vacuum to reproduce 
f n = 93 MeV, m n = 140 MeV as well as the s-wave ir-n scattering phase shift in the one-loop 
order: one of such parameter sets is —p 2 = (284MeV) 2 , A/47T = 5.81 and h 1 / 3 = 123MeV 
40| . This will be used through out this report. 

In Fig]6|, shown are the spectral functions Pn^) in the n (cr)-channel at finite T calculated 
in the 0(4) linear a model: Two characteristic features are the broadening of the pion peak 
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(Fig.3(A)) and the spectral concentration at the 2 m n threshold (a; ~ 2 m n ) in the a-channel 
(Fig.3(B)). The latter may be measured by the 2y spectrum from the hot plasma created in 
the relativistic heavy ion collisions H.0. 



FIG. 6. Spectral functions in the n channel (A) and in the a channel (B) for T=50, 120, and 
145 MeV with the parameter set shown in the text |4(|. p CT (u;) in (B) shows only a broad bump at 
low T (a), while the spectral concentration is developed as T increases, (a)—>(b)—>(c). 


B. Spectral enhancement in nuclear matter 


The idea of the spectral enhancement near the 2m n 

B 


the finite density system 


threshold can be also applied to 
A surprise was that such enhancement can be seen even 
near the nuclear matter density if the chiral condensate decreases by 30 % or so. To see 
this explicitly, let us consider again the 0(4) linear cr model to demonstrate the essential 
idea of the near-threshold enhancement. The Lagrangian adopted in |{43 
Gell-Mann-Levy model with a minor modification 


is the standard 


£mn — + V’W — 9 $)ip + 


( 20 ) 


where M = a + ir ■ tt, M 5 = cr + i'y §t ■ 7T, and ip is the nucleon field. In the mean-held level, 
medium effect originating from the nucleons can be incorporated by making the following 
replacement; 


(a)o -> (a) = $(p){cr)o, (21) 

where (-)o ((■)) denotes the vacuum (nuclear matter) expectation value, and p (po) is the 
baryon density (nuclear matter saturation density) .0 <&(p) may be parameterized as < h(p) = 
1 —Cp/po with C = 0.1 — 0.3 at low density. The resultant spectral function in the a channel 


2 p in the following corresponds to ns in previous sections. 
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is shown in Fig.[7j (left panel). The spectral enhancement just above the 2 m n threshold 
similar to the case at finite T (Fig.^j(B)) can be seen when (a) changes by 25 % from its 
vacuum value. This mean-field treatment was later improved by including the p-h and A -h 
contributions to the pion propagator f44[]. In this case, the spectral strength spreads into 
the energy region even below 2 r m 7T , but the qualitative feature of the enhancement at < lm v 
still remains as shown in Fig.|7] (right panel) where a = 1 — m*(p = po)/m*(p = 0). 



FIG. 7. Left panel : The spectral function in the u channel for several values of = (a) /(cio) 
with the same parameter set of the linear a model as Fig. (|6j) . This figure is taken from jH|. Right 
panel : The phenomenological spectral function in the a channel with the p-wave effects such 
as p-h and A -h excitations. The solid line corresponds to the vacuum case and the dashed lines 
correspond to the in-medium case at p = po with several values of a(= 1 — m* (p = po)/m*(p = 0)), 


One can also study the in-medium 7r-7r amplitude using the same model 
unitarized tc-tt amplitude in the I=J =0 channel reads 


45]. A simple 


T(s) 


1 

T t ^(s)-iQ(sY 


( 22 ) 


where T tree (s) is the in-medium tree-level amplitude with the mean-field replacement (pl|), 
and O(s) is the phase space factor defined as 0 


1 I 

0(s) = 0(s - - —1. (23) 

In Fig.|] (left panel), the in-medium tt-tt cross section cr^(s; p) is shown in the arbitrary 
unit, in which the same parameter sets as in Fig.|6j and Fig.[j] are used. Again a large 
enhancement of the cross section near the threshold is seen as the baryon density increases. 
In Fig.|] (right panel), shown is the cr T7r (s; p) relative to its vacuum value defined as 


3 Here we have changed the definition of T(s) by factor 2 from that in [45]. 
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( 24 ) 


p _ < ^7r7r('S, p) 

P 0) 

A large enhancement near 2m n threshold in Fig.|]can be understood as follows: As the 
partial restoration of chiral symmetry takes place and (a) decreases, m* approaches 2 m^, 
which implies that T^ e (s ~ 2m,) tends to be suppressed, ffence the s-dependence of the 
full inverse amplitude T _1 (s) just above the 2i r threshold is governed by the imaginary part 
@(s)/2, which causes the near-threshold enhancement of T(s). 




FIG. 8. Left panel: In-medium tt-tt cross section calculated in the linear a model in the 
mean-field approximation [45]. The same parameter set with Fig.(]6],^) are taken. Right 
panel : In-medium enhancement of the cross section relative to its vacuum value in the range 
2 m n < v's < 500 MeV @. 


At this point, a natural question to ask is whether the near-threshold enhancement 
obtained in the 0(4) linear a model arises also in the non-linear models or in the models 
which do not have a explicitly ||5]]. Furthermore, if it is the case, what kind of vertices 
in the non-linear chiral Lagrangian are responsible for the enhancement? To study these 
problems, it is best to start with the standard polar parameterization of the chiral held ||45|| , 
M = a + ir ■ tt = ((a) + S)U with U = exp {if ■ </>//*)• Here, (a) is the chiral condensate 
in nuclear matter as before, while /* is an “in-medium pion decay constant”. The nucleon 
held is also transformed from 0 to the chiral invariant held N. 

Now let us take a heavy-scalar ( S ) limit and heavy-baryon (N) limit simultaneously. 
These limits can be achieved by taking A,g —> oo with g/X and (<r)o = f n being hxed. 
In this limit, the heavy scalar held S may be integrated out and the following effective 
Lagrangian is obtained: 


h_ 

U 

where all the constants take their vacuum values: = (cr) 0 , m 2 a = A(<t)q/3 + m 2 , and 

ttin = g(cr)o. Note that gf 7X /2m 2 a in front of NN approaches to a finite value 3g/2Xf n in the 
heavy limit, thus it cannot be neglected. In (^5|), i s the standard p -wave n-N coupling 
and ■ • • denotes other higher dimensional operators which are not relevant for the present 
discussion. 


TMt/t + [/]]+ 44 + 


(25) 


jCoff — 


W - (Trpudun 


2 ml 
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In the uniform nuclear matter, NN in eq.(|25|) may be replaced by p. This leads to a 
reduction of the vacuum condensate; f n = (cr) 0 —> (a) = (a) 0 (1 — gp/ = f*. Then the 

proper normalization of the pion held in nuclear matter should be 0' = (<j)/ f n ) ■ /*. Thus, the 
origin of the near-threshold enhancement in the heavy limit can be ascribed to the following 
new vertex: 


£new — —Trry- NNTr[dUdtf]. 

7T 


(26) 


Because this vertex is proportional to the scalar-isoscalar density of the nucleon, it affects 
not only the pion propagator but also the interaction among pions in nuclear matter. In 
Fig I, 47T-IV-IV vertex generated by C new is shown as an example. Note that the vertex in 
Fig.|9] acts to enhance the Tin attraction in the I—J— 0 channel in nuclear matter. Despite 
its important role, this vertex has not been considered so far in the calculations of the Tin 
scattering amplitudes in nuclear matter in the non-linear approaches. 



FIG. 9. An example of the new Att-N-N vertex generated by £ nev 
and the dashed line represent the nucleon and the pion, respectively. 


The solid line with arrow 


Here it should be pointed out that the vertex eq.(p6|) has been known to be one of 
the next-to-leading order terms in the non-linear chiral Lagrangian in the heavy-baryon 
formalism [[hj. In fact, Thorsson and Wirzba [[L7|] have derived a general in-medium chiral 
Lagrangian by taking the mean-held approximation for the nucleon held; 


<r> = (f +1 p ) Tr[ara/t] 


+ (k- 


C_2 

2 

f 2 

' J 7T 


9a 


16m tv 

Cl 


) p Tiid 0 Ud„U'} 


(27) 


+ ( : j + YP)Tr(C/t x + x , f7). 


Note that - disappears in >177) in the mam-field level because of the isospin symmetry in 


"■kN 

nuclear matter, although it may contribute beyond the mean-held approximation as pointed 
out in |48|. As is shown in J|6|, the coefficients Cg 2,3 are related to the irN sigma term, the iso¬ 
spin even S'-wave nN scattering length, and the nucleon axial polarizability, respectively. 
We just quote here the numbers although they have large potential uncertainties: ci = 
—0.87GeV -1 , c 2 = 3.34G6V- 1 , and c 3 = -S^SGeV” 1 . 

By comparing eq. (P5|) with eq. (P7|) , we observe that the new vertex of the form 
NNTr[dUdU'] appear in both cases with the same sign. In eq.(^5|), we have only the Lorentz 
invariant terms because we have integrated out only the scalar meson S. On the other hand, 
eq.(|27|) contains more general 0(3) invariant terms because it potentially contains the effect 
of heavy excited baryons |^6| . It is of of great importance to make an extensive analysis of 
the mr interaction in nuclear matter with the new Att-N-N vertex. 
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C. Behavior of the a pole in 2nd Riemann sheet 


We emphasize here the relation of the near-threshold enhancement and the softening 
of the fluctuating mode in the a channel in nuclear matter. We shall show that such a 
fluctuating mode can be characterized by a complex pole of the unitarized amplitude T(s) 
in eq.(p2|). 

To see this in a simplest way, let us first take the heavy scalar limit ( m a —> oo) and the 
chiral limit (m T —> 0). In this case, the tt-tt correlation is generated only by the standard 
lowest order chiral amplitude in eq. (|22| ) with T tree (s) = s/f 2 . Then one readily finds a pole 
of T(s) in the lower half plane as [|45| 


\J Spole = V87T{cr) (1 - i), 

= (466 — 466i)y^— MeV. 
\°/o 


(28) 


This complex pole represents a broad damping-mode dynamically generated by the nn 
multiple scattering in the I=J =0 channel. Since (a) decreases in nuclear matter, the pole 
moves toward the origin in the complex y^-plane ||9| and thus causes an enhancement in 
the mi cross section as shown in Fig|| The existence of such softening in the complex plane 
and its connection with the near-threshold enhancement of the spectral function was first 


shown in |31|J33|] using the Nambu-Jona-Lasinio model. 

We shall now see that these qualitative features, i.e., the existence and the softening of a 
complex pole do not change in the more realistic case with the real part of the rescattering 
amplitude included properly. We start again with the non-linear model with m a —»■ oo. The 
T matrix resumed in the inverse amplitude method |5(J reads 


T(s) = T 2 /(l + T 2 g(s)), 


(29) 


where T 2 is the lowest order chiral amplitude (= s/f 2 in the chiral limit) and g(s) is the 
pion loop integral; 


g(s) =i j 


d A q 


(2vr) 4 q\ 


7nl + ieqi — ml + ie' 


(30) 


with q± — q ± p/2 and s = p 2 . We notice that Img(s) = —0(s), which is given in 
The real part of g(s) is logarithmically divergent, so one needs a renormalization or use the 
one-subtracted dispersion relation, which introduce an arbitrary constant. Actually, this 
program can be nicely formulated in the N/D method with the dimensional regularization 
15| . However, as is also shown in [HJ, the resultant formula can be obtained by a simple 


cut-off method provided that the cut-off value is properly chosen. The result in the chiral 
limit is given by [15 m 


9( z ) = 7^{a SL ~ Log(-/i 2 /^)}, 


167T 2 


for z in the first sheet, and 

9(z) = 


{a SL - Log (-pi 2 /z) + 2ni}, 


167r 2 


(31) 


(32) 
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for z in the second sheet (Im z < 0). Here Log z is the principal value of the logarithm; 
it has a cut along the real negative axis and — 7 r < Arg[Logz] < 7r. The constant a SL is a 
cut-off dependent constant and is determined so that the empirical phase shift in the scalar 
channel is reproduced |5T|; a SL ~ —.75 ± .2 with fi = 760 MeV. The pole position is given 


by the following dispersion equation 


t 2 ‘W = -sW- 


(33) 


The result is shown in Fig.5: When f n = 93 MeV, the pole exists at -y/s = 461 — 341z MeV. 
As f n =(a) is decreased, the pole moves toward the origin. It means that the mass and the 
width of the a mode decreases as the chiral symmetry is restored; thus we confirm that the 
a is the soft mode of the chiral transition. 



FIG. 10. The shift of the pole position of the T matrix in the chiral limit as / 7r =(o') decreases 
toward 0. 


V. POSSIBLE EXPERIMENTAL SIGNATURE 


A. Softening in the scalar channel 


The CHAOS collaboration measured the n + n ± invariant mass distribution ± in the 
reaction A(ji + , 7r + 7r ± )A / with the incident pion momentum p v + = 399 MeV and the targets 
A= 2 H, 12 C, 40 Ca, 208 Pb [|52f . They observed that the yield for M^ + ^_ near the 2m n threshold 
is close to zero for A = 2, but increases with increasing A. They identified that the 7r + 7r _ 
pairs in this range of M^ +7r _ is in the I=J =0 state. 

Although this experiment was originally motivated by the possible strong ti-tt correlation 
in nuclei due to many-body effects |)3) , state of the art calculations based on the conventional 
many-body theoretical approach without incorporating the effect of partial restoration of 
chiral symmetry do not reproduce the sufficient enhancement in the near-threshold region 
in the I=J =0 channel |7I|]. However, once the effect of partial chiral restoration in nuclei 


suggested in |13| is incorporated together with the conventional approach, the near threshold 
enhancement enough to explain the experimental data may be obtained |F|]. So far such an 
analysis is still in the qualitative level. Quantitative calculations based on both the linear 
and non-linear approach discussed in the previous section are called for. 
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For the experimental confirmation of the threshold enhancement seen in |52[, measure¬ 
ments of 27r° and 27 final states with hadron/photon beams off the heavy nuclear targets 
are necessary. Those channels are free from the p meson background inherent in the 7 r + 7 r _ 
measurement. One of such experiments measuring 27r° from the reaction A(n~, 7 r° 7 r°)kl / with 
p n - = 408 MeV and A— H, 2 H, C, Al, Cu has been done by the Crystal Ball Collaboration 
55fl . They report some qualitative difference of the n-n invariant mass distribution from the 
CHAOS data. For making a better comparison between the data from different detectors 
with different acceptance, the following “combined ratio” is useful ||56|| , 




a 


A 


(M wn )/a, 


a 


N 


(M nn )/cr, 


A 
T 
N ’ 


(34) 


where cr^ ( 07 ?) is the measured total cross section of the ( 7 r, 27 t) process in nuclei (nucleon): 
This ratio represents the net effect of nuclear matter on the interacting (7T7r)/ = j = o system. 
As shown in Fig.|TT|, it has recently been argued that both the CB data and the CHAOS data 


are consistent and gives an enhancement near the 2?7 i n threshold |57j, although the statistics 
in the CB data is relatively poor. For more detailed discussions, see the contributions to 
this Proceedings from CHAOS and CB groups. 
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FIG. 11. The composite ration as a function of the 7T-7T invariant mass for the CB data (filled 
diamonds) and CHAOS data (open diamonds). The target is 12 C in both cases. This figure is 
taken from 


We list some of the possible experiments related to the in-medium a in the following. 


• Measuring of 2y’s from the electro-magnetic decay of the a or (7r7r)/ = j =0 in nuclear 
matter, although the branching ratio is small, is an interesting alternatives to the 27T 
measurement because of the small final state interactions. In such an experiment, one 
needs to fight with large background of photons mainly coming from 7T°s. Nevertheless, 
if the enhancement is prominent, there is a chance to find the signal. 
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When a has a finite three momentum, one can also detect dileptons through the scalar- 
vector mixing in matter owing to the violation of the charge conjugation symmetry at 
finite density [0 


cr 


7 


e + e . (See also, 0.) 


The inverse process 7 + A —> 2tt + X can be also used to produce the a or 
system by the electro-magnetic probes. Such experiments with electro-magnetic probes 
have been planned and are being performed in SPring -8 [|59| . 


We remark that the nuclear reactions such as (d, 3 He) and (d, t) are also useful to 
explore the spectral enhancement and/or the formation of the cr-meson nuclei [j6C 


in the case of the deeply bound pionic nuclei and the possible production of 77 - 

0 


as 

or 


amnesic nuclei 


The incident kinetic energy E of the deuteron in the laboratory 


system is estimated to be l.lGeV < E < 10 GeV, to cover the spectral function in 
the range 2m 7r < uj < 750 MeV. 


B. Medium effect in other channels 


The dilepton spectrum measured in the heavy ion collisions such as Pb-Au collisions both 
at high energy (158 GeV/A) and at low energy (40 GeV/A) show a sizable enhancement of 
the e + e _ yield below the p-meson peak 0. This may or may not be related to the partial 
chiral restoration in nuclear medium originally proposed in 0] (for the recent review, see 
|64|]). I 11 Fig.[l2], the recent CERES/NA45 data for 40 GeV/A collision are shown. The E325 
experiment at KEK [0] measured e + e~ pairs from the p -A collision at 12 GeV. The similar 
enhancement over the known source and combinatorial background as CERES is seen in the 
mass range of about 200 MeV below the p-u peak for A = Cu (Fig.|T3|). 

The deeply bound pionic atom has proved to be a good probe of the properties of the 
hadronic interaction deep inside of heavy nuclei 0[. It is suggested that the anomalous 


energy shift of the pionic atoms (pionic nuclei) owing to the strong interaction could be 
attributed to the decrease of the effective pion decay constant f* (p) at finite density p which 
may imply that the chiral symmetry is partially restored deep inside the nuclei 
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FIG. 12. Normalized e + e _ invariant mass spectrum in the Pb-Au collision at 40 A GeV. Each 
contribution from the known hadronic sources are shown by thin solid lines, and their sum is shown 
by thick solid line. This figure is taken from the second reference in (62|. 


[events / 50MeV/c 2 ] [events / 50MeV/c 2 ] 




FIG. 13. e + e” invariant mass spectrum in the p-A reaction, (a) is for the C and polyethylene 
targets and (b) is for the Cu target. The solid lines in both figures are the result of the best-fit 
including known dilepton sources and the combinatorial background. Excess over the solid line is 
seen below the p-u peak in the case (b). 


VI. SUMMARY AND CONCLUDING REMARKS 

This report may be summarized as follows. 

The QCD phase diagram should be considered in mi-T-ns space where rrii (i = u,d,s ) 
being the current quark mass, T ( ns ) is the temperature (baryon density), the recent 
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extensive studies based on the lattice QCD simulations and effective theories of QCD indeed 
shows it quite probable that the chiral phase transition takes place for T c ~ 170MeV and 
71b ~ several x no, although the precise order of the transition in the real world is not clear 
yet. 

The order parameter of the chiral symmetry breaking is ( qq ). Furthermore, its amplitude 
and phase fluctuations are the urcr and the pion, respectively, which play a central role in 
hadronic world. They mix strongly through the process a which may account for 

the fact that the physical a is elusive in the real world. Irrespective of the introduction 
of the ur-cr, it has become more and more plausible in recent years that the strong n-n 
attraction produces a light a pole in the 2 nd sheet of the tt-tt scattering amplitude in the 
I=J =0 channel. Constraints from chiral symmetry, analyticity and crossing symmetry are 
key ingredients for the confirmation of the light cr-pole. Such a scalar-isoscalar fluctuation 
has also phenomenological importance in modern hadron physics. 

If T and/or ns of the system increases, the intimate connection between the QCD vacuum 
structure and the fluctuation modes becomes even more apparent. In fact, the light a meson 
acts as a soft mode and the complex cr-pole in the second sheet has characteristic shift toward 
the real axis in association with the partial restoration of chiral symmetry (qq) —> 0 . 

In cold nuclear matter, a sizable decrease of the quark condensate, as large as 30 % 
at normal nuclear matter density, may take place according to the low-density expansion 
of (qq). If this is the case, probing the scalar-isoscalar fluctuation in cold nuclear matter 
using heavy nuclei is promising for exploring the partial restoration of chiral symmetry. The 
softening of the complex cr-pole leads to a strong enhancement of the spectral function and 
the 7T-7T cross section in the I=J =0 channel near the 2 m n threshold. This phenomena can 
be explicitly demonstrated e.g. in the linear and nonlinear realization of chiral symmetry 
provided that the possible reduction of the quark condensate or is taken into account. 

The possible relevance of the above theoretical considerations and the ( 7 r + 7 r _ )/ = j = o data 
(the CHAOS collaboration) and the ( 7 r° 7 r°) data (the Crystal Ball collaboration) is suggestive 
and should be studied further. Independent tests using the electro-magnetic probes as well as 
the mesic-nuclei will be also useful to unravel the realization of chiral symmetry in hadronic 
matter. 

This work is partially supported by Grants-in-Aid of the Japanese Ministry of Education 
Science and Culture (No. 12640263 and 12640296). 
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